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Abstract: We show how such important features of QCD as chiral symmetry breaking 
or the formation of a mass-gap can be directly traced from QCD sum rules for two point 
functions assuming, in the large number of colors limit, exact duality between the operator 
product expansion and the spectrum described by linearly (or nearly linear) rising Regge 
trajectories as predicted by string theory. We see how the presence of chiral symmetry 
breaking is intimately related to confinement in this scenario, as expected from general 
arguments, and how Regge trajectories change when chiral symmetry is broken. As a 
result the whole meson mass spectrum can be parametrized with a good accuracy by the 
constant f n only, thus realizing the program proposed by Migdal some time ago. 
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1. Introduction 

In the present work we consider QCD sum rules for the vector (V), axial- vector (A), scalar 
(S), and pseudoscalar (P) channels in the large- N c limit, the planar limit of QCD [1] and 
match the corresponding operator product expansions (OPE) with the predictions of the 
QCD effective string theory. This approach has been used by several authors previously, 
(see e.g. [2-7]). In a sense, this work may be partly viewed as a simplified version of the 
analysis performed in [2] but with a rather different accent. 

First of all, we show how the parameters of the effective string mass spectrum can be 
extracted from the sum rules for the two-point functions in a very simple way and without 
involved numerical fits (this is the main distinction with [2]). Second, we try to give a 
clear physical interpretation of all results appearing in the analysis. We see how in this 
framework chiral symmetry breaking and confinement (i.e. the creation of a mass gap) 
are intimately related. Finally we provide in this rather well defined theoretical setting an 
implementation of Migdal's program whereby all physical quantities in theories like QCD 
are expressed, modulo simple numerical factors, in terms of a basic scale. f n is this scale 
in the present case. 

The paper is organized as follows. In section 2 we introduce the relevant formulas. 
Sections 3 and 4 are devoted to the analysis of linear string-like mass spectrum for vector 
and scalar channels correspondingly. In sections 5 and 6 this analysis is extended to the 
non-linear spectrum. The results obtained are interpreted in section 7. Section 8 is devoted 
to a comparison with previous results. Final section 9 contains some concluding remarks. 
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2. General formalism 



Let us briefly remind the reader of the general formalism for deriving QCD sum rules in 
the large- N c limit. In this limit the two-point correlation functions of quark currrents are 
saturated by an infinite set of narrow meson states with the quantum numbers of these 
currents, i.e. they can be represented in Euclidean space as follows 

U J (Q 2 ) = [ d A xe^(qTq(x)qr q (0)) planar = £ J"?^]. r + D J Q + D(Q\ (2.1) 
J n H + m j\ n ) 

J = S,P,V,A; T = ^75,7^,7^75; D , D x = const. (2.2) 

For simplicity we do not write the Pauli matrices but the SI// (2) symmetry is understood. 
The last two terms are required for renormalization of infinite sums. On the other hand, 
their high-energy asymptotics are given by the OPE [8-10] (we consider the chiral limit) 

11 W ) = To - ? 7T2 ^ To 7TZ I n 7rQ!s— TTfi - ' ( 2 - 3 ) 



12vr 2 Q 2 12vr Q 4 9 Q 
U^(Q 2 ) — -^-Q 2 In— + — • - M! (2 4 ) 



where 



?> P = -7, e A ' 5 = 11, (2.5) 



and we have defined 



n^(Q 2 ) = (SfivQ 2 + QpQv) U V ' A (Q 2 ). (2.6) 



The symbols (G 2 ) and (gg) denote the gluon and quark condensate, respectively. The 
residues are parametrized as follows 

Z VA (n) = 2F 2 :A (n), Z s , P {n) = 2G% P {n)m\ P {n) . (2.7) 

with Fv : A(n) being electromagnetic decay constants and Gp{n) are related with the cor- 
responding weak decay constants (see [2] for details) . 

The sum rules simply follow from comparison at each power of Q~ 2 of the OPE (2.3) 
and (2.4) with the sum in Eq. (2.1) after summing up over resonances (in a chiral invariant 
way) and subtracting infinite constants which are irrelevant for our purposes. 

3. Sum rules for linear spectrum: Vector case 

Phenomenology tells us that the linear mass spectrum predicted by effective string theories 
(such as the ones represented by the Veneziano amplitude, the supersymmetric string, or 
the Lovelace-Shapiro amplitude) is a very good approximation to the real world [11]. In the 
present section we shall adhere to this type of spectrum, namely assume that the spectrum 
is strictly linear and see what are the implications of this assumption. Thus consider the 
ansatz 

m 2 (n)=m 2 + an, F 2 (n) = const, n = 0,1,2,... (3.1) 
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Here the universal slope a is proportional to the string tension T: a = 2ttT. After renormal- 
ization (subtraction of infinite constant) and taking into account the 7r-meson contribution 
(or more generally, the possibility of a zero-mass state in the pseudoscalar and axial-vector 
channels) we have 



n y,A^_2/^ ,y 2F h _2f 2F^ A ( C? + m* \ 

n (Q ) - -^Saj + X, q 2 + m 2 + on " qT^j {—^— ) > ( 3 - 2 ) 



where we have used the Kronecker symbol 5 a b to indicate that the hypothetical Goldstone 
boson contributes only to the axial-vector channel due to PCAC. 
The ^-function has an asymptotic representation at z S> 1 

k=i 



there i?2fc denote the Bernoulli numbers. 
Let us introduce the notation 



2 

m 



x= . (3.4) 



a 



Performing the procedure outlined above one arrives at the following sum rules 

N c F] 



24vr 2 o 

p2 



(3.5) 



= f^Aj ~ Fj (x - 1/2) , (3.6) 
^(G 2 ) = Fja(x 2 -x + l/6), (3.7) 

^7ra s (^) 2 = -\F]a 2 x (x - 1/2) (x - 1) . (3.8) 

From Eq. (3.5) it follows that the quantities Fj are directly expressed through the 
slope a. Let us substitute them into Eq. (3.8). The result is 

cJ (g a \2 

-— 7TQ s 167r 2 ^- = x(x - 1/2) (x - 1) . (3.9) 

N c a d 

Now we note that the quantity (qq) 2 /a 3 w (220/1100) 6 ~ 10 -4 is very small in the real 
world. Thus the theory has a natural small parameter. More precisely, since the quantities 
C J /3 and ira s are of order one, the small parameter is 

16^ 2 M^ w 0.01. (3.10) 

The existence of this small parameter in the sum rules seems to be tightly related to that 
of the chiral perturbation theory, i.e. the ratio squared of m n to the masses of light mesons 
in other channels. A kind of "perturbation theory" can be developed here. In the zero 
order the sum rules become extremely simple providing nevertheless a good approximation 
to the real world. In the rest of the paper we will consider only this case. 
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In this approximation Eq. (3.9) takes the form 



x (x - 1/2) (x - 1) ~ 0, (3.11) 

and only two dimensionful parameters appear: The slope a which parametrizes the strength 
of the gluon interaction and, hence, must be directly related to the gluon condensate {G 2 ) 
(indeed we see that from (3.7)), and the 7r-meson weak decay constant f w parametrizing 
chiral symmetry breaking 1 . 

Let us investigate the solutions of the above set of equations. First of all, equations 
(3.5)-(3.7) and (3.11) have one trivial solution (even without any approximations) corre- 
sponding to the case where the theory is weakly coupled (non-confining), without bound 
states and condensates: Fj = a = but the ratio in Eq. (3.5) is non-zero and give the 
parton-model result. All non-perturbative effects are absent in this case. 

Let us now analyze the possibility of a nontrivial 2 chirally symmetric solution with 
f n = (here of course (qq)=0 and Eq. (3.11) is strictly true). From Eq. (3.6) we have only 
one possibility, namely x = 1/2, which is consistent with Eq. (3.11). Hence the chirally 
symmetric linear spectrum would be 

m 2 (n) = ^+an. (3.12) 

Substituting the solution in Eq. (3.7) and making use of Eq. (3.5) we obtain a relation 
between the slope and the gluon condensate 

numerically resulting in ^-(G 2 ) ~ — (370 MeV) 4 that is in agreement with the standard 
value in QCD, ^{G 2 ) = (360 ± 20 MeV) 4 , in absolute value, but not in the sign, which is 
an obvious inconsistency. 

In fact the gluon condensate is proportional to the QCD vacuum energy e vac (see 
discussions in section 7). As QCD is an asymptotic free theory the vacuum energy turns 
out to be negative for the stable vacuum [13]. We draw attention to the quite interesting 
fact that the solution x = 1/2 is nothing but the exact minimum of the function 3 e vac ~ 
— (G 2 ) ~ x 2 — x + 1/6. Hence, in the vector channels the chirally symmetric case for the 
linear spectrum corresponds to the theory near the maximum of this function (since the 
sign is opposite to that of vacuum energy). 

Consequently, the chirally symmetric solution we have just been exploring cannot 
exist for the linear spectrum if confinement holds. Thus, the property of confinement in 
the large- N c limit automatically results in Chiral Symmetry Breaking (CSB), i.e. f n ^ 0. 
This agrees with the Coleman- Witten theorem [14]: if confinement persists in the planar 
limit of QCD, chiral symmetry is inevitably broken. As we will see later, the analysis of 



1 Throughout the paper we use the value of / T in the chiral limit [12], f„ = 87 MeV. 
2 I.e. with a mass gap. 

3 The vacuum energy (or gluon condensate) in our approach turns out to be the function of intercept, 
which is quantized, i.e. this energy is also quantized. 
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the chirally broken case necessarily requires the introduction of non-linear corrections to 
the spectrum in any case. 



Staying within the framework of the linear spectrum, we note that the chirally symmet- 
ric spectrum (3.12) for the case of p- mesons agrees with the one derived from the so-called 
Lovelace- Shapiro (LS) amplitude [15] for the reaction tt + it — > it + n, which had a great 
phenomenological success in the 60's 



r(l - a p (s) - a p {t)) 



with ctp(s) = \ + | . This amplitude predicts equal masses for the ground vector and 
scalar mesons, i.e. the solution x = 1/2 would extend (if one were to believe the LS 
amplitude) to both cases. An important property of this amplitude is the presence of an 
'Adler zero' [16], i.e. the fact that A(s,t) — > as s,t — > 0. This is actually a consequence 
of Goldstone theorem so there is some logical uneasiness here since we have assumed so 
far that chiral symmetry is unbroken, but being the presence of the Adler zero a necessary 
but not sufficient condition for CSB one can find a scape route from contradiction. 

The generalization of this LS amplitude to reactions of the form ir+A — > B+C was per- 
formed in [17]. The generalized LS amplitude describes the axial- vector and pseudoscalar 
states. It turns out that the amplitude derived in [17] predicts x = 1 for axial-vector 
and pseudoscalar mesons. This is numerically inconsistent with the sum rules (even if we 
include f n ) as they stand and, most importantly, with chiral symmetry restoration at high 
energies (see e.g. [2]). The conclusion is that LS, in spite of its undeniable phenomeno 
logical success, is unable to reproduce the sum rules for two point functions. This rules it 
out as a possible model for QCD even leaving aside the fact that it does not derive from 
any known string theory. 

Coming back to Eqs. (3.5)-(3.8) and having ruled out the possibility of a chirally 
symmetric solution to the sum rules we have to look for solutions with f w / 0. We know 
from current algebra that this implies a non-zero value for the quark condensate. However, 
as we have discussed previously, the numerical value for the physically relevant quantity 
167T 2 (qq) 2 1 a 3 is very small, so we can still use the approximate equation x(x — \)(x — 1) ~ 0. 
If x = \ is ruled out because it corresponds to the chirally symmetric solution, can we 
perhaps consider the other two solutions x = 0orx = l? The first one is obviously ruled 
out due to sign inconsistencies. The second one is viable, but inconsistent with the sum 
rules unless we accept deviations from strictly linear trajectories. As we do not know 
how CSB distorts the linear spectrum, we will introduce these effects phenomenologically 
through non-linear corrections to the string-like spectrum. Obviously chiral symmetry 
restoration at high energies requires adopting the same value of x in all channels [2,6] and 
we will consider non-linear corrections in sections 5 and 6. 
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4. Sum rules for linear spectrum: Scalar case 



It is convenient to make the rearrangement [2]: 
G 2 (n)m 2 (n) 



n(Q z 



Q2 _|_ m 2( n ) 



+ D + DiQ z 



^2G 2 (n)+D c 



Q 2 



E 



2G 2 (n) 
Q 2 + m 2 (n) 



-D 1 



(4.1) 



where the infinite constant Dq and D\ subtract the infinities coming from the sums over 
resonances. After renormalization we have 



n 



S,P 



Q 2 + m 2 + an 

n=0 



^ 6PJ + 



2G 2 SP Q 2 



Q 2 + m 2 



(4.2) 



Let us consider again the chirally symmetric case. The pion residue vanishes in this case 
and it will be taken into account in section 6 where the chirally broken case is considered. 
Performing the standard procedure (just as in the vector channels using the same notation 
x) one obtains the sum rules 



16vr 2 



G 2 ; 



-G 2 ja (x 2 



2£" 



x + 1/6) , 
- V2) (x - 



I) 



(4.3) 
(4.4) 

(4.5) 



-7ra s {qq) 2 = -G 2 a 2 x(x 

The general requirement for the spectrum is the chiral symmetry restoration at high 
energies (i.e. equal x for both channels). There are three solutions with these properties: 
x = 0, 1/2, 1. Substituting these solutions in Eq. (4.4) and making use of Eq. (4.3) we 
obtain the relation between the slope and the gluon condensate 



a. 



7T 



(G 2 ) 



a. 



(G 2 



N c a 2 

24vr 2 ' 

N r a 2 



for x 



1 

2' 



for x = , 1. 



(4.6) 



(4.7) 



vr ' ' 12vr 2 ' 

The first case is in agreement with the spectrum of the LS amplitude for the scalar channel 
(but not for the pseudoscalar channel) and the gluon condensate is reproduced with the 
correct sign this time (see discussion after Eq. (3.13)). The lightest pseudoscalar has the 
mass, m 2 = a/2, in agreement with the well known fact that if no spontaneous CSB occurs, 
then pion stays massive. The appearance of the other two solutions is a consequence of 
the fact that an analogue of Eq. (3.6) is absent in the scalar channels. The vacuum energy 
is evidently larger on these solutions, hence, they are unphysical. It is interesting to note 
that the gapless solution that appears is among the energetically unfavorable ones. The 
existence of this solution expresses another exact result: If a rigorously chirally symmetric 
theory possesses a massless pseudoscalar bound state then there is a degenerate scalar 
meson partner (see, e.g., [18]). The fact however that the value of (G 2 ) disagrees in sign 
with Eq. (3.13) is nevertheless another sign of the inconsistency of the chirally unbroken 
solution. 
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5. Sum rules for non-linear spectrum: Vector case 



According to the discussions at the end of section 3 we are forced to introduce a non-linear 
spectrum. We will derive the minimal non-linear realization of the spectrum consistent 
with the sum rules. For that we have to make some hypotheses about what represents the 
leading contribution in the non-linear case. First of all we know nothing from experiment 
about the residues of radially excited mesons. Numerical calculations show that they differ 
negligibly from the linear case [2]. Thus for simplicity we do not introduce non-linear 
corrections for the residues of excited states and for ai-meson as well (below we show that 
F ai nicely agrees with experiment without any corrections). So, except for the p- meson, 
we denote 

Fl = F\ = F 2 . (5.1) 

We note also that the mass of p-meson is in the real world substantially below 1 GeV, i.e. 
the typical CSB scale. Consequently it is natural to expect that this meson is affected by 
CSB to a larger extent than other vector mesons. We shall make the hypothesis that CSB 
induces a shift in the residue of the p meson (and in no other vector or axial mesons) . We 
shall later attempt to justify this hypothesis. Then 

F 2 = F 2 + F 2 . (5.2) 

Now let us introduce the non-linear mass spectrum in the form 

m 2 j(n) = m 2 + an + 5j(n), (5-3) 

where the corrections always can be represented as follows: 

Sj(n) = aAjfj(n), /j(0) = 1. (5.4) 

Here the dimensionless constants Aj are supposed to be due to CSB, i.e. they are propor- 
tional to f 2 /a, and fj(n) are decreasing functions of n only They have to vanish at least 
exponentially in n according to [2], or polynomially according to [6] (see discussion in sec- 
tion 8), but in this work we will not adhere to any particular model for these functions. If 
the non-linear corrections provide a rather big contribution to the masses of ground states 
but relatively small one to those of excited states (as we expect) then the parametrization 
(5.4) automatically will not be very sensitive to a concrete choice of fj(n). 

An important assumption of the present analysis is that the slope a has no corrections. 
This is a universal quantity which does not feel CSB. Any corrections to the spectrum are 
encoded in the value of intercept m 2 and in the non-linear contribution 5{n). Such a picture 
is well justified phenomenologically. 

The contribution of the non-linear mass corrections to the correlator is (for the time 
being we denote 5 = 5j (n) and omit the general factor) 



1 1 



1 + 



Q2 _|_ m 2 _|_ an _|_ J Q2 _|_ m 2 _|_ an \y Q2 _|_ m 2 _|_ an ^ 
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The /9-meson additional contribution is 



Q 2 + m 2 + aA v Q$ \Q 
Let us introduce the constants 

CO OO CO 

C J ^/j(n), = TTnfj(n), C 2 J EE^/j(n). (5.7) 

n=0 n=0 n=0 

After summation over n and making use of notations (5.1), (5.2), and (5.7) sum rules 
(3.5)-(3.8) become 



24vr 2 ~ a 



N c F l 



(5.8) 



f 2 8 AJ + F 2 5 V j = F 2 (x — 1/2) , (5.9) 
^T a {G 2 ) = x 2 - x + 1/6 - 2Af , (5.10) 

= ar(x-l/2)(x-l)-3A^, (5.11) 
where the corrections A/ are non-zero due to CSB and they are given by 

F 2 

A{ = AjCi + j^(x + A v )5 V j, (5.12) 
A 2 J = Aj {2Cix + 2C( + CiAj) + ^(x + A V ) 2 5 VJ . (5.13) 
From Eq. (5.9) we get 

F = U. (5.14) 

Thus from the same equation one has 

1 f 2 1 24vr 2 /^ 
^=2 + ^ = 2 + ^T- (5 - 15) 

At this stage we need some additional input to advance. For reliability this input should 
be a quantity which is experimentally known with a good accuracy. Except f w (which is 
already an input) the best candidate for such a quantity is F p = 154 ± 8 MeV. We note 
now that this number is very close to y/3 f n ~ 151 MeV. Thus, within our accuracy, we 
may put F p = y/3 as input. Making use of Eqs. (5.1) and (5.14) one obtains 

F = "v/2/tt ~ 123 MeV, (5.16) 

which is in agreement with the experimental value for the decay constant of ai-meson 
(see (5.1)): F ai = 123 ± 25 MeV. From Eq. (5.15) we immediately get the solution for the 
intercept, x = 1, and the relation 

a = f/fl (5-17) 

in very good agreement with phenomenology. This formula happens to be the result re- 
ported in [19] based on completely different grounds. It is curious to note that rela- 
tion (5.17) can be obtained from combining the slope of the LS amplitude, a = 2m 2 (we 
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remind that the slope has no corrections in the present analysis), with a known formula for 
the mass of /j-meson (5.18) (both relations were also found in [5]). In principle, Eq. (5.17) 
could be taken as input from the very beginning, all other results are then automatically 
reproduced. This way would be natural in the sense that considering the CSB case we 
should insert a spectrum of string with CSB incorporated. Such strings were considered, 
e.g., in [20]. Further development of this approach resulted in a qualitative derivation 
of Eq. (5.17) in [19], where the scales of the Goldstone boson physics and of the string 
dynamics are supposedly related through Eq. (5.17) by the chiral anomaly. 

It is thus quite remarkable that CSB forces a jump from x = \ to x ~ 1. A physical 
interpretation of the intercept is that it is a typical mass of ground state, hence, it should 
be tightly related to the scale of CSB, Acsb — 4^/^. Accidentally or not, they are identical, 
m 2 ~ a ~ A^ SB . A deeper understanding of this coincidence might lead to an independent 
derivation of Eq. (5.17). 

From the last sum rule Eq. (5.11) one can see that at the solution x = 1 the contribu- 
tion due to the linear part of the spectrum (the first term in the r.h.s.) is exactly cancelled. 
As a result the corrections turn out to be directly related to the (unwritten) quark con- 
densate in the l.h.s. This is in a full agreement with our initial assumption concerning the 
nature of the corrections. 

The slope a has a purely gluonic nature and, consequently, its value should not change 
dramatically after CSB. Thus, one may substitute Eq. (5.17) into Eq. (3.12) and obtain 
the mass of ground state in the chirally symmetric case: 

m 2 j(0) = fffl (5.18) 

This is nothing but a formula for the mass of p-meson obtained some time ago within the 
Borel QCD sum rules [21] and rederived in some other types of sum rules (see, e.g., [5, 
22]). This formula also usually holds when matching different models to the vector meson 
dominance and appears all over the literature (see, e.g., [23,24]). Relation (5.18) reproduces 
the experimental mass of p-meson with a unexpectedly high accuracy (we remind that 
f w = 87 MeV in the chiral limit). Thus, we should conclude that the mass of the ground 
vector meson practically does not shift after CSB. Then we immediately obtain the constant 
Ay parametrizing corrections in the vector channel 

24tt 2 f 2 1 

Av = = -I. (5.19) 

N c a 2 

Substituting Eqs. (5.14), (5.19), and (4.6) (which can be rewritten as ^-(G 2 ) = F 2 a) into 
Eq. (5.10) we obtain the estimate 

< = (5-20) 
From the same equation for the axial-vector channel we have the relation 

A A C£ = 1. (5.21) 

For a rough estimate of the ai-meson mass we may accept a kind of universality, Cq ~ Cq , 
which leads to m 2 ^ « 1.1a « (1150 MeV) 2 . Exactly the same value was obtained in the 
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framework of QCD sum rules [9]. It seems that a better estimate is hardly possible within 
the sum rules, let alone the fact that in the present work we have accepted the large- 
N c limit. This example shows that the mass corrections are indeed relatively small for 
the mesons with the mass above the CSB scale. Even for the ai-meson they constitute 
approximately 10% which competes with the accuracy of the large-iV c counting. This is 
what we expected from the very beginning. Keeping in mind that the constants Cf are of 
order 1 thanks to the assumed fast convergence of functions fj(n), we could independently 
recover in this way the approximate solution x = 1 in the CSB phase from Eq. (5.11) for 
the axial-vector case (knowing from Eq. (5.9) that x = and x = 1/2 are not solutions). 

6. Sum rules for non-linear spectrum: Scalar case 

As in the previous section we will be interested in the minimal non-linear realization of the 
spectrum consistent with the sum rules. Thus, we do not introduce non-linear corrections 
for the residues of excited states. Let us denote 

Gl = & P = Gl (6.1) 

For the ground states parametrization (2.7) should be changed 

Zj(0) — Zj(0) + Zj. (6.2) 

Here Zp is the 7r-meson residue which cannot be taken into account within parametriza- 
tion (2.7) and Zg reflects a contribution to the residue of ground scalar state. Introducing 
non-linear corrections to the masses in the form (5.4) and repeating the simple calculations 
outlined in the previous section one arrives at the sum rules 



N c _G 2 Q 



16vr 2 



(6.3) 



^(G 2 ) = -G 2 Q a {x 2 -x + 1/6) + Zj, (6.4) 



Zs&js 



= x (x - 1/2) (x - 1) - 3(z + Aj) \AjC J q + ^ . (6.5) 

Motivated by the previous results, we assume that the linear part of the spectrum does 
not induce CSB. The physical solutions of the system are then x = 1 for the scalar channel 
and x = 1 or x = for the pseudoscalar one. In both cases we obtain the relation 

Using the current algebra result, 

Z P = 2% (6.7) 

J TV 

and phenomenological values for the gluon condensate and the slope one can estimate the 
quark condensate 

(GO) ~ -(170 MeV) 3 . (6.8) 
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Taking into account our approximations, this estimation is selfconsistent. Combined Eqs. (6.6) 
and (6.7) pressupose a certain relation between the vacuum condensates and hadronic pa- 
rameters, which is similar to a result in [25]. For x = 1 one has Ap = —1 and the n-th 
pseudoscalar meson is a chiral partner for the n-th scalar meson. In the second case Ap = 
and the (n + l)-th pseudoscalar meson is a chiral partner for the n-th scalar meson. In this 
case there is no non-linear correction to the pseudoscalar spectrum and the 7r-meson has 
no chiral partner. 

For the scalar channel we have from Eqs. (6.5) and (6.6) 

A S C* = -1. (6.9) 

In order to estimate the mass of the ground state we can, like in the axial-vector channel, 
assume an approximate universality for the non- linear contributions, Cq ~ Cq . Making 
use of Eq. (5.20) this leads to m|(0) pa 0.7a sa (920 MeV) 2 . The relative contribution to 
the residue of ground state due to non-linear corrections is of order Zs/Zs(0) = 0.1 0.2, 
i.e. much less than for the p-meson. 



7. Discussions 



Let us give an intuitive interpretation which lies behind the results obtained. First of all, 
the gluon condensate is a measure of vacuum energy £ V ac- To be more precise we remind 
that the vacuum energy and the nonperturbative part (n.p.) of v.e.v. (G 2 ) are related by 
the anomaly in the trace of energy-momentum tensor [26] 

4e vac = (ep n . p . = ^(G 2 ) n .p. + 0(a) + ■ ■ ■ • (7.1) 

The term O(a) is the contribution of quark polarization effects. In [27] the effective po- 
tential for 0^ was constructed in the tree approximation. It was shown that the minimum 
of this potential (vacuum energy) must be negative, i.e. as the /^-function is negative only 
a positive gluon condensate provides a stable vacuum. 

In our analysis it is instructive to regard mesons as the spectrum of excitations given 
by the effective normalized quark vacuum "potential energy" U, where 

U = - — — 77j-^7> (7.2) 

7T a/(47T 2 ) 

and the gluon condensate (G 2 ) depends on the vacuum energy e vac . 

Let us take the vector and axial-vector channels as an example and discuss the three 
solutions of the sum rules mentioned in the paper. The trivial solution corresponds to the 
weakly-coupled QCD where the three- and four-point gluon interactions are suppressed 
and we effectively have U m i n = like in the QED. Since the vacuum is perturbative there 
are neither condensates nor resonances. 

The chirally-symmetric (CS) solution corresponds to fluctuations over the perturbative 
vacuum. The potential energy U has a local maximum in this point, U max = m 2 , which 
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results in a negative gluon condensate. As we still calculate the excitations over the zero- 
level, there appear the mass gap Acs = m2 p = a/2. 

After the CSB the potential energy recovers its minimal value, U m i n = — m 2 This 
new minimum turns out to be lower than the perturbative one by the value 8ir 2 f 2 which is 
equal to m 2 As a result the mass-gap enlarges by two times, Acsb = 2m 2 = a = (4irf w ) 2 . 
It is exactly equal to the generally accepted value for the gap squared of spontaneous chiral 
symmetry breaking A CSB ! 

Let us try to understand the " special" nature of p- meson in the framework of presented 
approach. In the perturbative vacuum we always have photons which are the only massless 
colorless particles in the Standard Model. However, in the nonperturbative vacuum this 
vector particle "feels" the "deep" below {7 = through its virtual quark-antiquark com- 
ponents. As a result the quark and antiquark in the virtual pair can acquire a dynamical 
(constituent) mass, converting the photon into a massive vector particle. In other words, 
we can observe the photon over this deep as a vector particle with the mass m p . Since the 
height of hump over zero-level and that of deep below this level are equal (the perturba- 
tive vacuum is distorted symmetrically in the both directions), the p-meson mass is not 
shifted after the CSB. The described effect in a natural way results in the enhancement 
(compared to other mesons) of probability for creation of this meson from vacuum, i.e. the 
positive shift of p-meson residue. This interpretation provides an intuitive picture about 
the mechanism of vector meson dominance. 

In the axial-vector channel the situation is completely different. The ai-meson pole 
is shifted in the CSB phase supposedly due to the ir — a± mixing, i.e. the ai-meson 
acquires additional mass due to a very intensive interaction with the Goldstone bosons. 
The mechanism seems to be the following. In the chirally symmetric phase both masses 
are equal, m\ x = m 2 = a/2. In the CSB phase these particles become maximally mixed. 
If two particles are mixed, rotation of their fields to a physical basis is known to result in 
repulsion of the initial masses. If the mixing is maximal the repulsion is then also maximal 
and in our case leads to m 2 ai = a/2 + a/2 = a, m 2 = a/2 — a/2 = 0. This interpretation 
is in a full agreement with old results of current algebra [28,29]. Namely, if the chiral 
symmetry were unbroken the ground states would fall into irreducible representations of 
SU(2) (g) SU(2) chiral group, with the masses being equal insight a given representation. 
After CSB these representations get mixed. In particular, the following relation holds: 
m 2 cos 2 ip + m 2 ai sin 2 ip = to 2 , where ijj is a mixing angle between pseudoscalar and axial- 
vector representations. This angle was determined from the p width to be around 45°, i.e. 
the mixing is maximal. In the chiral limit the relation converts then into the Weinberg 
formula. It is interesting to note that the scalar and vector representations practically do 
not mix after CSB. It is another manifestation of the fact that p-meson mass does not shift 
after CSB. So it does for the ground scalar meson. In practice, however, this state could 
have glueball, four-quark and strange components which hamper the prediction of its mass 
on the base of considered sum rules. 

Thus the underlying reason for the appearance of mass shifts in the CSB phase is the 
formation of constituent quarks which interact with the Goldstone bosons in a different 
way in channels with different quantum numbers. This effect is drastic for the ground 
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states, the higher is the radial excitation the less important is the CSB phenomenon and 
the chiral symmetry gets restored [6,30]. 

Let us emphasize that the formulas like m p = v^^/tt make sense only in the CSB 
phase. This point sometimes causes a confusion in the literature (see e.g. [22]). The 
mentioned relation by no means signifies that in-medium mass of p-meson changes in 
response to changing f n , it can not be an order parameter for phase transitions like /„-. 
The primordial formula is m p = \J a/2, i.e. this mass has a gluonic origin. What we should 
expect at high temperatures or baryon densities is a decreasing m ai up to the value m p , the 
latter being unchanged in the first approximation. At the same time should increase 
up to the same value m p . Finally, all masses of ground states become equal to \Jaj2 in 
the chirally restored phase. 

8. Relation to other work 

Concerning the nature of non-linear corrections to the string-like spectrum we would like to 
mention an interesting paper by Weinberg [29]. In that approach any mass matrix can be 
written as the sum of chirally invariant term and chirally non-invariant one. The latter term 
appears due to a non-trivial matrix of axial couplings and destroys the chiral symmetry 
in the spectrum. Our analysis suggests a one-to-one correspondence between the chirally 
invariant linear part of the spectrum and the first term, and the non-linear contribution 
and the second term in [29]. For the higher excitations the matrix of axial couplings is 
expected to gradually get trivial [31] resulting in suppression of the second term. This 
scenario realizes an expectation (see, e.g., discussions in [32] and references therein) that 
restoration of chiral symmetry in the hadron spectrum is related to decoupling of Goldstone 
bosons from highly excited states. 

In some models based on sum rules at large- N c [4, 5, 7] the intercepts in spectrum of 
chiral partners are different. This discrepancy with our results (and those of [2,6]) can be 
confusing. We note, however, that such type of spectra are usually obtained using some 
additional assumptions. In [7] both models were designed with the aim to obtain a non- 
zero dimension- two gluon condensate. The price to pay is, e.g., that the well-motivated 
relation (5.17) does not hold. In the model in [4] the expression for low-energy constant Lg 
was imposed as "another one" sum rule and then fine-tuned by a special prescription for the 
cut-off in the infinite sums. As a result this constant comes from perturbative logarithm, 
i.e., it is identically zero for any finite number of resonances (lack of a smooth transition 
from finite number of states to infinite one). The ideas in [5] are somewhat close to those 
of present paper, but in that paper there is no natural transition to the chirally symmetric 
case and the chosen way to introduce the non-linear corrections is rather questionable. 

In the literature exist two different estimations for the rate of chiral sym- 
metry restoration. In [6] the power-law minimal fall-off was obtained from 
matching to OPE without any convergence requirements for the sum rules, 
while in [2] the fall-off was estimated to be exponential after imposing the con- 
vergence for the generalized Weinberg sum rules. At present it is difficult to 
say which variant is realized in nature. 
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Finally we would like to mention a possible relation of our work to the so-called 
AdS/QCD models which are believed to provide a precise correspondence between con- 
formal gauge theories (like QCD in the ultraviolet region) and some superstring theories 
(see, e.g., [33] and references therein). These models can naturally incorporate chiral prop- 
erties of QCD, in particular, the Gell-Mann-Oakes-Renner formula was derived in [34] and 
relation to some effective models with hidden gauge bosons was demonstrated in [24]. As 
was noticed in [6] and developed in [35] they are often dual to the Migdal approach to 
the QCD sum rules [36] based on Pade approximation of the two-point correlators. Thus, 
QCD sum rules in the large- N c limit and AdS/QCD models seem to be complementary 
approaches. In this respect an interesting result was reported in [37]. Namely, at certain 
constraints on the infrared behavior of AdS/QCD one can obtain the linear spectrum for 
the vector and axial mesons with the same slope. The next successful step in this direc- 
tion would be a derivation of AdS/QCD model which reproduces equal intercepts to be 
consistent with chiral symmetry restoration at high energies. Another challenge for the 
AdS/QCD approach would be derivation of the relation for the string slope, Eq. (5.17). 

9. Conclusions 

In the late 70's before the publication of Shifman-Vainstein-Zakharov paper [8] A. Migdal 
formulated the program of calculation of the whole spectrum for mesons made of light 
quarks as a function of some unique dimensional parameter \i: M; L = fiCi, where the factor 
Cj depends on the discrete input parameters of QCD such as the number of colours and 
flavors [36]. This formula was supposed to be possible due to small values of bare quark 
masses which can be neglected (except the Goldstone bosons, of course). 

The approach presented in this paper is, in a sense, a concrete realization of Migdal 
program. Indeed, the spectrum depends on the constant /„-, the slope of trajectories, 
and the gluon condensate. These three parameters are related by Eqs. (4.6) and (5.17). 
Consequently, we may choose only one of them to parametrize the whole spectrum. The 
most natural choice is evidently f n . The first correction to the spectrum comes from taking 
into account the quark condensate squared just as was supposed by Migdal in the same 
paper! The extension of this analysis to mesons with other quantum numbers and to the 
SUf(3) group is a rather technical work. 

Acknowledgments 

We are grateful to A. A. Andrianov for fruitful discussions and enlightening comments. 
The work was supported by CYT FPA, grant 2004-04582-C02-01, and by CIRIT GC, 
grant 2001SGR-00065. The work of S.A. was also supported by RFBR, grant 05-02-17477. 
S.A. thanks the Ministry of Education and Science of Spain for financial support. 

References 

[1] G. 't Hooft, Nucl. Phys. B 72 (1974) 461; 
E. Witten, Nucl. Phys. B 160 (1979) 57. 



-14- 



[2] S.S. Afonin, A. A. Andrianov, V.A. Andrianov and D. Espriu, J. High Energy Phys. 04 
(2004) 039; hep-ph/0403268. 

[3] A. Bramon, E. Etim and M. Greco, Phys. Lett. B 41 (1972) 609; 
M. Greco, Nucl. Phys. B 63 (1973) 398; 
J.J. Sakurai,P/M/s. Lett. B 46 (1973) 207; 

N.V. Krasnikov and A. A. Pivovarov, Yad. Fiz. 35 (1982) 1270; Phys. Lett. B 112 (1982) 
397; 

A. L. Kataev, N.V. Krasnikov and A. A. Pivovarov, Phys. Lett. B 123 (1983) 93; 
S.G. Gorishnii, A.L. Kataev and S.A. Larin, Phys. Lett. B 135 (1984) 457; 
E.A. Tainov, Z. Physik C 10 (1981) 87; 

B. V. Geshkcnbcin, Yad. Fiz. 42 (1985) 991; 

B.V. Geshkenbein, Sov. J. Nucl. Phys. 49 (1989) 705; 

M. Shifman, At the Frontier of Particle Physics/Handbook of QCD, (ed. M. Shifman), World 
Scientific, 2001; hep-ph/0009131; 

M. Goltcrman and S. Peris, J. High Energy Phys. 01 (2001) 028; 
S.R. Beane, Phys. Rev. D 64 (2001) 116010; 
S.R. Beane, Phys. Lett. B 521 (2001) 47; 

Yu.A. Simonov, Phys. Atom. Nucl. 65 (2002) 135; hep-ph/0109081; 

M. Golterman, S. Peris, B. Phily and E. de Rafael, J. High Energy Phys. 01 (2002) 024; 

A. A. Andrianov, V.A. Andrianov and S.S. Afonin, hep-ph/0212171; 

M. Goltcrman and S. Peris, Phys. Rev. D 67 (2003) 096001; 

V.I. Shevchenko and Yu.A. Simonov, Phys. Rev. D 70 (2004) 074012; hep-ph/0406276; 

A.L. Kataev, Phys. Atom. Nucl. 68 (2005) 567; hep-ph/0406305; 

O. Cata, M. Golterman and S. Peris, J. High Energy Phys. 05 (2005) 076; 

J.J. Sanz-Cillero, Nucl. Phys. B 732 (2006) 136; hep-ph/0507186; 

V.D. Orlovsky, hep-ph/0510192; 

S.S. Afonin, A. A. Andrianov, V.A. Andrianov and D. Espriu, Int. J. Mod. Phys. A 21 (2006) 

885; hep-ph/0509144; 

S.S. Afonin, hep-ph/0606291. 

[4] O. Cata, M. Golterman and S. Peris, hep-ph/0602194. 

[5] S.S. Afonin, Phys. Lett. B 576 (2003) 122. 

[6] M. Shifman, hep-ph/0507246. 

[7] E.R. Arriola and W. Broniowski, hep-ph/0603263. 

[8] M.A. Shifman, A.L Vainstein and V.I. Zakharov, Nucl. Phys. B 147 (1979) 385, 448. 

[9] L.J. Reinders, S. Yazaki and H.R. Rubinstein, Nucl. Phys. B 196 (1982) 125. 

[10] M. Jamin and M. Miinz, Z. Physik C 60 (1993) 569; hep-ph/9208201. 

[11] A.V. Anisovich, V.V. Anisovich and A.V. Sarantsev, Phys. Rev. D 62 (2000) 051502; 
V.V. Anisovich, Usp. Fiz. Nauk. 47 (2004) 49; hep-ph/0208123; 
D.V. Bugg, Phys. Rept. 397 (2004) 257. 

[12] J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 

[13] M.A. Shifman, A.L Vainstein and VI. Zakharov, Sov. Phys. JETP Lett. 27 (1978) 55; 
A. A. Andrianov, V.A. Andrianov, V.Yu. Novozhilov and Yu.V. Novozhilov, Phys. Lett. B 
186 (1987) 401; 

A. A. Andrianov and V.A. Andrianov, Z. Physik C 55 (1992) 435. 



- 15 - 



[14] S. Coleman and E. Witten, Phys. Rev. Lett. 45 (1980) 100. 

[15] C. Lovelace, Phys. Lett. B 28 (1968) 264; 
J.A. Shapiro, Phys. Rev. 179 (1969) 1345. 

[16] S. Adlcr, Phys. Rev. 137 (1965) B1022. 

[17] M. Ademollo, G. Vencziano and S. Weinberg, Phys. Rev. Lett. 22 (1969) 83. 

[18] A. H611, A. Krassnigg and CD. Roberts, Phys. Rev. C 70 (2004) 042203. 

[19] A. A. Andrianov, D. Espriu and A. Prats, AIP Conf. Proc. 756 (2005) 302; 
hep-ph/0412380. 

[20] J. Alfaro, A. A. Andrianov, L. Balart and D. Espriu, Int. J. Mod. Phys. A 18 (2003) 2501. 

[21] M.A. Shifman, A.I. Vainstcin and V.I. Zakharov, Phys. Rev. Lett. 42 (1979) 297. 

[22] E. Marco and W. Weise, Phys. Lett. B 482 (2000) 87. 

[23] E.R. Arriola and W. Broniowski, Phys. Rev. D 67 (2003) 074021; 

E. Megias, E.R. Arriola, L.L. Salcedo and W. Broniowski, Phys. Rev. D 70 (2004) 034031. 

[24] D. T. Son and M. A. Stcphanov, Phys. Rev. D 69 (2004) 065020. 

[25] S.S. Afonin, hep-ph/0603176. 

[26] J.S. Collins, L. Duncan and S.D. Joglekar, Phys. Rev. D 16 (1977) 438. 

[27] A. A. Migdal and M.A. Shifman, Phys. Lett. B 114 (1982) 445; 
M.A. Shifman, Phys. Rept. 209 (1991) 341. 

[28] F. J. Gilman and H. Harari, Phys. Rev. 165 (1968) 1803. 

[29] S. Weinberg, Phys. Rev. 177 (1969) 2604. 

[30] A. Le Yaouanc, L. Oliver, S. Ono, O. Pcne and J.-C. Raynal, Phys. Rev. D 31 (1985) 137; 
A. A. Andrianov and V.A. Andrianov, hep-ph/9705364; hep-ph/9911383; 
A. A. Andrianov, D. Espriu and R. Tarrach, Nucl. Phys. B 533 (1998) 429; 
A. A. Andrianov and D. Espriu, J. High Energy Phys. 10 (1999) 022; 
L.Ya. Glozman, Phys. Lett. B 475 (2000) 329; 

A. A. Andrianov, V.A. Andrianov and S.S. Afonin, hep-ph/0101245; hep-ph/0209260; 

V.A. Andrianov and S.S. Afonin, Yad. Fiz. 65 (2002) 1913; hep-ph/0109026; 

V.A. Andrianov and S.S. Afonin, Eur. Phys. J. A 17 (2003) 111; 

V.A. Andrianov and S.S. Afonin, J. Math. Sc. 125 (2005) 99; hep-ph/0304140; 

E.S. Swanson, Phys. Lett. B 582 (2004) 167; hep-ph/0309296; 

T. DeGrand, Phys. Rev. D 69 (2004) 074024; 

L.Ya. Glozman, A.V. Nefediev and J.E. Ribciro, Phys. Rev. D 72 (2005) 094002; 
L.Ya. Glozman, hep-ph/0411281; 

S.S. Afonin, Phys. Lett. B 639 (2006) 258; hep-ph/0603166; 
S.S. Afonin, hep-ph/0606310. 

[31] S.S. Afonin, hep-ph/0605102. 

[32] L.Ya. Glozman and A.V. Nefediev, hep-ph/0603025. 
[33] I.R. Klebanov, hep-ph/0509087. 

[34] J. Erlich, E. Katz, D.T. Son and M.A. Stephanov, Phys. Rev. Lett. 95 (2005) 261602. 



- 16 - 



[35] J. Erlich, G.D. Kribs and I. Low, Phys. Rev. D 73 (2006) 096001. 
[36] A. A. Migdal, Ann. Phys. (NY) 110 (1978) 46. 

[37] A. Karch, E. Katz, D.T. Son and M.A. Stephanov, hep-ph/0602229. 



-17- 



